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Abstract 

A finite difference computer model is developed to simulate the exposure statistics of a radio 
frequency buoyant antenna as it is towed in a three-dimensional random seaway. The model 
allows the user to prescribe antenna properties (length, diameter, density, etc.), sea conditions 
(significant wave height, development of sea), tow angle, and tow speed. The model then 
simulates the antenna-sea interaction for the desired duration to collect statistics relating to 
antenna performance. The model provides design engineers with a tool to predict antenna 
performance trends, and to conduct design tradeoff studies. The floating antenna envisioned 
is for use by a submarine operating at modest speed and depth. 
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Chapter 1 



Introduction 

1.1 Motivation 

As the United States Navy launches the last of the highly capable Seawolf class submarines and 
begins construction on the even more capable Virginia class submarines, it continues to find itself 
operationally limited by the technological barriers present with high data rate communications. 
The capability for high speed data transfer has escalated at an astonishing rate over the last 
few years; however, the submarine force faces unique challenges in fully utilizing high data 
rate systems which operate at higher frequencies. High frequency (UHF, SHF, or EHF band) 
radio signals are unable to penetrate below the surface of the water; so, while the submarine 
is submerged, radio communications systems are limited to low data rate receive-only modes 
transmitted at very low frequencies (ELF, or VLF band). The current mode of operation calls 
for the submarine to routinely come to periscope depth to raise a radio mast above the surface of 
the water, send and receive the necessary message traffic, then submerge to continue its mission. 
This required visit to periscope depth not only causes the submarine to suspend execution 
of its mission, but also reduces the overall stealth of the mission. Submarine operations 
are typically driven by the need for stealth, and stealth is best preserved by remaining fully 
submerged. Clearly, the communications goal of the submarine force is to be able to utilize 
the high frequency and high data rate communication circuits while retaining the stealth of 
operation at speed and depth. 



9 



1.2 Background 



This goal of high data rate coiiiinuiiication at speed and depth is the motivation of a U.S. Navy 
Advanced Technology Demonstration (ATD). Initially funded jointly by the Defense Advanced 
Research Projects Agency (DARPA) and the Office of Naval Research (ONR) in FY 96-97 as 
a feasibility study, the Buoyant Cable Array Antenna (BCAA) research became a DARPA and 
ONR Demonstration in FY 98-99. The BCAA prototype satisfactorily demonstrated receive- 
only capability with a field test in the Summer of 1999. The program w<is promoted to the level 
of Navy ATD in FY 00 and renamed the Multi-element Buoyant Cable Antenna (MBCA). The 
MBCA is a 3 year joint program which is a combined Navy ATD-funded and DARPA-fimded 
effort. This ATD, one of only 6 ATD’s funded by the Navy for fiscal year 2000, reflects the 
high priority the Navy places on submarine connectivity. The goals of the MBCA program 
include: 

♦ demonstrate 24 kbps 2- way communication at UHF over the Fleet Satellite (FLTSAT) 
communications system while the submarine is at a speed of 6 knots and at a depth of 
300 feet 

♦ demonstrate 1.5 Mbps data rates over UHF Line-of-Site 

• demonstrate L-band voice communications. 

The participants in the program include the Naval Undersea Warfare Center, Johns Hopkins 
University Applied Physics Laboratory, Sippican Inc., and MIT’s Lincoln Laboratory. Lincoln 
Laboratory is a federally funded laboratory administered by the Massachusetts Institute of 
Technology. Lincoln Laboratory, whose specialties include radar and electrical engineering, 
was chosen to overcome the technological challenge of designing a floating linear array radio 
frequency communication antenna module and the dissociated inboard electronics. However, to 
gain insight into the behavior of the antenna module in a random seaway, Lincoln Laboratory 
added Professor Jerome Milgram of the MIT Ocean Engineering department to the antenna 
development effort. Professor Milgram joined a working group that was established to estimate 
statistics on 

• the average fraction of aiitemia elements exposed over an extended period of time 
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• the temporal variation of the exposure statistics, including duration of reception dropout 
times 

• the spatial distribution of those exposed elements. 

This thesis work was undertaken to support the statistics working group by improving upon 
the first hydrodynamic model which was developed by Professor Milgram and Gary Ulrich. 
[9] [10] The improved model retains all of the capabilities of the original model and improves 
upon the original model by upgrading the wave field from a unidirectional wave field with the 
choice of only two tow directions, 0 or 180 degrees, to a multidirectional wave field with the 
choice of any tow direction from 0 to 180 degrees. This provides all angles due to port and 
starboard symmetry. 

1.3 Antenna Model 

The method used to determine the antenna behavior and the exposure statistics was to develop 
long term simulations from which exposure statistics could be obtained vice getting the exposure 
statistics from the theory of a non-linear process. The model, named PROTEUS^ performs 
the simulation by constructing a wave field and then modeling the antenna in the wave field 
and analyzing the resultant behavior. The development of the wave field model is described in 
Chapter 2 and the development of the antenna model is described in Chapter 3. The results 
of the improved model are presented in Chapter 4 and a detailed user’s guide is included as 
Appendix A. 
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Chapter 2 



Wave Field Model Development 

2.1 Linear, Random, Deep Water Wave Theory 

To simplify the model and to reduce computation time, the simulated wave field was assumed to 
be described by linear, random, deep-water wave theory. The linear theory aspect of the waves 
implies that the velocity potential and Huid velocity are constant from the mean free-siirface to 
the free-surface level. The linear theory assumptions represent a first order approximation in 
satisfying the free surface kinematic and dynamic boundary conditions. A complete derivation 
and discussion of linear wave theory can be found in Faltinsen and Newman. [4] [6] “Deep 
water” is defined to exist when A < 2d, where A is the wavelength and d is the water depth. 
The random aspect of the waves implies that 

• the instantaneous wave elevation from the mean free-surface is gaussian distributed with 
zero mean and variance 

• the waves are an ergodic process, which means that the wave statistics obtained from a 
single location are tlie same as the wave statistics olitained from an ensemble of locations 

• the waves are a stationary process, which means that the wave statistics are invariant 
with respect to time [7] [8] 

Because the hist condition is so restrictive, we relax it to allow the statistics to vary slowly 
with time, otherwise known as a weakly-stationary process. A summary of the deep water 
linear theory wave relations is provided in Table 2.1. 
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Quantity 


Formula 


Wave Number 
Circular Frequency 


II II 


Dispersion Relation 


g 


Velocity Potential 


(f) = sin{kx' — ut) 


Surface Displacement 


h = acos{ka/ — ut) 


Radian Frequency 


II 


Phase Velocity 


1'^ 

II 

o 


Group Velocity 


II 


Horizontal Fluid Velocity 


u = cos(fcx' — cot) 


Vertical Fluid Velocity 


lu = sin(A:x' — ut) 


Wave Energy Density 


E = \pga? 


Note: x' refers to a stationary 


X coordinate system. 



Table 2.1: Linear Theory Deep Water Wave Relations 



Assuming a lower wave frequency bound of 0.3 rad/sec^ which is reasonable for the vast 
majority of energy density spectra one would encounter, this translates into a minimum water 
depth of 342 meters. This implies that for depths shallower than about 350 meters, the 
deep water theory begins to break down, and a finite depth treatment of the waves must be 
considered. A complete derivation and discussion of finite depth wave theory can be found in 
Faltinsen and Newman. [4][G] 

Upon inspection of the linear theory formulas in Table 2.1, the waves are found to be 
sinusoidal. However, an observation of the sea surface shows that the waves display a quite 
complicated and random structure. This observed behavior can be numerically approximated 
by adding together a series of sinusoids with each discrete frequency component having the 
appropriate amplitude and a random phase. As the number of discrete frequency components 
increases, the numerical approximation approaches the continuous distribution of the actual 
sea surface. This discrete method of describing a two dimensional random seaway is presented 
in Equation 2.1 



N 

h{x\t) = '^^ai{uj)cos{ki{uj)x^ — uJit + (f)^) (2.1) 

i=l 



where 
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Bretschneider Energy Density Spectrum 




Figure 2-1: Typical Wave Energy Density Spectrum 

h{x^ y t) is the wave elevation, referenced to the mean free-surface 

(ii{uj) is the appropriate wave amplitude of each discrete frequency component 

ki{uj) is the wave number of each discrete frequency component 

Ui is the set of discrete circular frequencies contained in the modeled energy density spectrum 
(j)^ is the random phase assigned to each discrete frequency component 
and the term has been dropped, since evaluation occurs at the sea surface {z = 0). 
Using the equations found in Table 2.1, one can determine most of the terms on the right 
hand side of Equation 2.1: ki is determined from the dispersion relation, and is a random 
number with a uniform distribution on 0 to 27T. The determination of the appropriate wave 
amplitude, however, requires some discussion and a deeper understanding of wave theory. 

To determine the wave amplitude corresponding to a given wave frequency, one begins by 
defining the spectral energy density function, S{u)^ which represents the energy of the waves 
having frequency u, A typical energy density spectrum is shown in Figure 2-1. The units of 
S{u)) are defined to be in ni^-sec, so that when S{lo) is integrated with respect to cj, the result 
is the average energy of the random waves with respect to time. 
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Bretschneider Energy Density Spectrum 
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Figure 2-2: Determination of Component Amplitudes ai (uj) 



The spectrum shown in Figure 2-1 is said to be a single-sided spectrum, which refers to 
the fact that the frequencies are constrained to be positive. For a double-sided spectrum, both 
the positive and the negative frequencies are allowed, and the spectrum is defined such that 
S{cj) = S{—cj), In addition, the double-sided spectrum is half the magnitude of the single-sided 
spectrum; since, for constant sea severity, the area under the spectral energy density curve must 
be constant. The amplitudes, ai{u)^ of the discrete frequency components in Equation 2.1 are 
obtained from Equation 2.2 



where 

ai{uj) is the amplitude of the discrete band, Ao;, of waves with center frequency lj. 
Equation 2.2 is simply a statement that the square of the wave amplitude is a measure of 
the average energy of the waves, because both sides of the equation represent energy in units 




( 2 . 2 ) 



of m?. The method of obtaining the amplitudes Oj(a;) is shown in Figure 2-2. 
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Until now, S{uj) has been presented in terms of wave energy present at a given wave fre- 
quency. It can be shown from the Wiener-Khintchine theorem that S{lj) is related to the wave 
amplitude record at a fixed location in the ocean. [7] We define a wave record at a given 
location as Then the autocorrelation of x\t) is given by 

1 7 ^ 

R{t) = lim — / x’{t) x'{t + r) dt (2.3) 

i — ^oo ZJ. J 

and by the Wiener-Khintchine theorem, the spectral energy density function, S{uj), and the 
autocorrelation function, i?(r), form a Fourier transform pair. 



S{u) = 


1 


noo 

/ R{t) e-^'^dr 


(2.4) 




7T j 


' — OO 




1 


roo 




R{r) = 


1 S{uj) e-^^duj 

— OO 


(2.5) 



This method proves to be invaluable as it has the ability to convert observed wave records into 
spectral energy density functions. 



2.2 Energy Density Spectra Contained in the Model 

Numerous energy density spectral forms have been proposed by researchers attempting to char- 
acterize a particular sea wave environment. [4] [5] [7] [8] The spectral forms may characterize a 
particular geographic wave environment or they may characterize a particular storm environ- 
ment, but in all cases, they share a very similar form. The model incorporates the spectral 
formulation as a subroutine; therefore, it is completely arbitrary from the model viewpoint how 
the spectrum is defined. In order to provide some variability in the allowable spectral forms 
and to allow for troubleshooting and debugging of the model, four wave spectral energy forms 
were included in the model. Two of the spectral forms are actual derived models that are well 
accepted in the ocean design field. Though the two spectral models may appear very different 
at first glance, most spectral formulations share the basic form of Equation 2.6 [9] 

S{u;) = (2.6) 
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where 



B are constants to be deterniined 

Hs is the significant wave height, which is defined to be the average wave height of the 
one third highest waves. Sea severity is most commonly given in terms of significant wave 
height among researchers; however, significant data exists which gives sea severity in terms of 
a value called sea state. Appendix B provides a correlation amongst significant wave height, 
modal wave frequency, wind speed, and sea state. The two actual derived spectral formulations 
chosen for the model are the Bretschneider two parameter spectrum and the Ochi six parameter 
spectrum, given in Equations 2.7 and 2.8, respectively. [7] 



The other two models are the square wave spectrum and the single sine wave spectrum, 
given in Equations 2.9 and 2.10, respectively. These simple geometric form models were used 




(2.7) 




( 2 . 8 ) 



for troubleshooting and debugging of the model and to verify the results of the model. 



/• 



0.0 for 6 <0.01 



5 ( 6 ) = 



5.0 for 0.01 <h< 1.01 



(2.9) 



0.0 for 5>1.01 



5 (6) = 5(5 (6 - hs) 



( 2 . 10 ) 



where 
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Bretschneider Energy Density Spectrum 
(De\«lopment = 1.0 for all cases) 




Figure 2-3: Bretschneider Spectrum for Varying Sea Severity 




and hg is the spatial frequency of the sine wave. 



( 2 . 11 ) 



2.2.1 Bretschneider Two Parameter Spectrum 

The two parameters of the Bretschneider spectrum, Equation 2.7, are significant wave lieight, 
Hg, and modal frequency, ujm> Modal frequency is the frequency at which S{uj) reaches its 
maximum value, and significant wave height is as previously defined. For the special case of 
a fully developed sea, a relationship exists between and Hg, A fully developed sea, also 
referred to cis a fully arisen sea, is one in which the wind imparts energy to the waves at a 
rate equal to that at winch viscous damping dissipates their energy and wave- wave interactions 
transmit energy from one frequency to another. Designating this modal frequency as 
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Bretschneider Energy Density Spectrum 
(Hg = 3.0 m for all cases) 




Figure 2-4: Bretschneider Spectrum for Varying Degrees of Development 



have 



^mpA 




( 2 . 12 ) 



which establishes the energy equilibrium characteristic of a fully arisen sea. The two require- 
ments for a fully developed sea are 



• sufficient distance over which the wind is blowing over the water, or fetch 

• sufficient length of time the wind has been blowing, or duration. 

If the fetch or duration are insufficient, the modal frequency will be greater than that given 
by Equation 2.12; if the winds are subsiding, the modal frequency will be less than that given 
by Equation 2.12. This is best understood by considering the mechanism by which waves are 
generated. The wind causes ripples, or capillary waves, on the sea surface, which is a short 
wavelength, high frequency phenomenon. Thus as a storm is building, wave energy input occurs 
at high frequency, and the energy gradually is transferred through wave-wave interactions to 
lower frequencies. Likewise, as a storm is subsiding, the high frequency waves attenuate most 
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rapidly due to viscous effects, and the result is a skewing of the spectral energy toward lower 
frequencies. By introducing a term we will call the development factor, Dev^ we can quantify 
the degree of sea arousal as given in Equation 2.13. [9] 

~ Dev • (2.13) 

Using this convention, values of Dev <1.0 correspond to a decaying sea, and values of Dev >1.0 
correspond to a building sea. Figures 2-3 and 2-4 show representative plots of S{u;) as significant 
wave height and development, respectively, are varied. Note that in Figure 2-4, that while 
changing the development shifts the location of peak spectral energy density, it keeps the area 
under the curve,and thus the total energy, constant. 

2.2.2 Ochi Six Parameter Spectrum 

The Ochi six parameter spectrum is a generalization of the Bretschneider two parameter spec- 
trum in two respects: 

• two spectral peaks are present, each with their own modal frequency and significant wave 
height, and their contributions add together 

• a shape factor. A, can be applied to each spectral peak to skew the curve right or left, 
independent of the significant wave height or modal frequency. [9] 

In the form given in Equation 2.8, each parameter is an independent variable. However, 
by analyzing over 800 wave spectra from the North Atlantic, Ochi obtained a relationship for 
the six parameters in terms of the measured significant wave height, Hg- [7] The result of this 
analysis is a family of spectra having a 95% confidence; however, the results given in Table 

2.2 represent a single “most probable” spectrum [7]. The model implements this single most 
probable form for obtaining the six parameters of the Ochi spectrum. A sample plot of the 
Ochi most probable spectrum, showing its characteristic bimodal shape, is shown in Figure 2-5. 
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Parameter 


Value 


Hsi 


0.84//^ 


Hs2 


0.54//^ 


^m\ 




^T712 


1.1 


Ai 


3.00 


A 2 


1.54e“°°'^2//5 



Table 2.2: Ochi Six Parameter Spectrum: Most Probable Values 



Ochi Energy Density Spectrum 
(Single Most Probable Spectrum for = 3.0 m) 




Figure 2-5: Ochi Six Parameter Single Most Probable Spectrum 
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Square Wa\e Energy Density Spectrum 




Spatial Frequency (m'^) 



Figure 2-6: Square Wave Spectrum 

2.2.3 Square Wave Spectrum 

The square wave spectrum was added to the model to provide for verification of the resultant 
multidirectional wave field spectral analysis. By providing a more easily recognizable input 
spectrum, the goal was to generate a more easily recognizable output spectrum. A plot of the 
square wave spectrum is shown in Figure 2-6. 

2.2.4 Single Sine Wave Spectrum 

The single sine wave spectrum was added to the model to provide for an efficient recognizable 
input spectrum which could be used for troubleshooting the frequency related components of 
the multidirectional wave field. By providing a more easily recognizable input spectrum, the 
goal was to generate a more easily recognizable output spectrum. A plot of the single sine 
wave spectrTim is shown in Figure 2-7. 
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Single Sine Wave Energy Density Spectrum 
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Figure 2-7: Single Sine Wave Spectrum 

2.3 Unidirectional Modeling Methods 

As mentioned earlier, the wave field in the previous model consisted of a unidirectional random 
seaway with the choice of only two tow directions. The wave field in the current model makes 
use of the unidirectional method to construct the multidirectional wave field, so it is necessary 
to completely understand the unidirectional method developed by Ulrich before moving on 
to the development of the multidirectional method. The discussion and derivation of the 
unidirectional method that follows comes from Ulrich. [9] 

The sea domain variables of interest for the unidirectional wave field are 

• sea elevation, h 

• fluid horizontal velocity, u 

• fluid vertical velocity, w. 

Two methods were employed in arriving at the sea domain properties of the unidirectional 
wave field, each with certain advantages and disadvantages. The first method will be called 
the Sine Wave Superposition (SWS) method, and the second method will be called the Fast 
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Fourier Tiansform (FFT) method. The source code for both of these computer programs can 
be found in reference [10]. 



2.3.1 Sine Wave Superposition (SWS) Method 

As previously mentioned in Equation 2.1, the sea elevation can be obtained by a superposition 
of sine waves having a random phase and an amplitude prescribed by Equation 2.2. In like 
manner, the fluid velocities can be calculated by the superposition method. The formulas for 
/i, u and w in an earth-fixed reference frame are provided in Equations. 2.14 - 2.16. 



N 



h{x\ t) 


= ^ ai{u) cos{ki{u) X - Uit -f 

N 


(2.14) 


u{x\t) 


= ai{uj) uj cos{ki{u) x^ — uit -f (p^) 

i=\ 

N 


(2.15) 


w{x' ^t) 


“ CJ sin{ki{uj) x^ — uit -f (p^) 

2=1 


(2.16) 



By invoking a coordinate transformation and the dispersion relation from Table 2.1, one arrives 
at the forms appropriate to the coordinate system moving with the antenna. These are presented 
in Equations 2.17 - 2.19. 



N 



h{x^t) = y^ ai{uj) cos 
i=l 
N 

u{x, t) = Oi{uj) U) cos 
i=l 
N 

w{x,t) = a; sin 



ki{u) X - [iJi )t + (pi 

9 



i=l 



ki{(jj) X — [Ui jr + (p. 



ki{uj) X — {uji )t + <f>^ 

9 



(2.17) 

(2.18) 
(2.19) 



The spectral energy density function S{u) used in this formulation was a single sided spectrum. 
The advantages of the SWS method are 



• simplicity in coding fluid property subroutines 

• good resolution for S{u) at lower frequencies, where most of the wave energy resides 
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Time- Frequency 
Domain 


Space-Spatial Fre- 
quency Domain 


f=l/T 


b= 1/L 


U = 27 t / 


k = 2nb 


k = {2nf)'^lg 


U = y/27Tgb 



Table 2.3: Comparison of Fourier Transform Domain Relationships 



• frequency meshing (Ao;) is unrelated to the spatial meshing (A.x), allowing both antenna 
element length and wave energy frequency resolution to.be arbitrarily chosen. 



2.3.2 Fast Fourier Transform (FFT) Method 



This method takes advantage of the fact that Equations 2.14 - 2.16 are a Fourier series, and 
thus the functions /i, w,and iv are the inverse Fourier transform of some expression. Since the 
calculations to determine the sea domain properties occur at a fixed point in time, the transform 
variables are not u and i, rather they become b and x. Table 2.3 summarizes the analogous 
relationships between the time-frequency domains and the space-spatial frequency domains. 
Because the FFT algorithm used for the model assumed the presence of both positive and 
negative frequency components, it was necessary to use a two sided spectrum for this method. 
In order to retain the same amount of energy (i.e. area under the spectral energy density curve) 
as the transition from one-sided to two-sided spectra is made, it is necessary that the formula 
for the component amplitudes be adjusted from that given in Equation 2.2. In addition, it is 
convenient to incorporate the random phase, 0^, into a complex component amplitude, q. With 
these changes in mind, the formulas providing the sea domain properties are given in Equations 
2.20 - 2.24. 



h{x) 

u{x) 

w{x) 



N 



Re |^Ci(6i)exp i{2itbiX -f {2nbiU - y/2ngbi) t) | 

Re I ^ Ci{bi) \/^2ngbi exp |^z( 27 r 6 ia: -h {2nbiU - ^^271 gbi) i) j 



Re \/ 2ngbi exp |^2(27r62X + {2nbiU — y/2ngbi) ^)j 



( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 
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Re {cj(6i)} = \I-S{bi) Ab cos(</>J 



(2.23) 



Im{ci(6i)} = \l-S{bi)Ab s\n{(f>i) 



Then the Fourier transform pairs become 



where 



and 



is defined to mean 



i=l 



H{b) h{x) 



N 



H{b) = ^ h{x) exp {~i27rbx) 

1 ^ 

— 7T~ {i^T^bx) 

/.7T ^ 






By setting the complex component amplitudes, Cj, to obey the relation 



(2.24) 



m 


O h{x) 


(2.25) 


\/2TTgbiH{b) 


u{x) 


(2.26) 


-iy/2ngbiH{b) 


w(x) 


(2.27) 


N 

11(b) = J2ci{b) exp 


27rbiU — >/27r^6j^j 


(2.28) 



(2.29) 



(2.30) 

(2.31) 






(2.32) 



that is, the positive frequency components are the complex conjugates of the negative frequency 
components, then the imaginary parts of Equations 2.20-2.22 are zero, and the Re{ } can be 
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dropped from the equations. Finally, tlie wave field is constructed by performing an inverse 
FFT on the one-dimeiisional complex arrays, represented by the left side of Equations 2.25 - 
2.27, which contain the Fourier coefficients describing h,u, and xv. 

When converting between the spectral energy density function as a function of u to another 
variable, such as 6, the basic relations of Equations 2.33 and 2.34 are useful. 

S{uj)(Lj - S{b)db (2.33) 

S{b) - 5(0.)^ (2.34) 

An important feature of the FFT method is the linkage of the spatial frequency domain 
meshing (A6), and the sea spatial meshing (Ax). The key to the FFT method is to make the 
antenna spatial meshing, As, match the sea domain spatial meshing. Ax. This is accomplished 
by using the relations given in Equations 2.35 - 2.38, and ensuring that the sea Ax equals the 
desired antenna As. 



Ab = 


1 . 

Y — ^min 

^sea 




(2.35) 


bmax ^ 


AtgeaAb ^ sea 

2 2Ls^a 


1 

~ 2 Ax 


(2.36) 


^min ~ 


V^27T^6niin 




(2.37) 


^max “ 


gbmax 




(2.38) 



The primary advantage of the FFT method is computational speed. The FFT method 
executes approximately 10 times faster than the SWS version. A further increase in execution 
speed is accomplished by not computing the FFT every time step. At. Rather, the sea data 
(h,u, ix) is computed at nAt intervals, and linearly interpolated in between. The integer n is 
determined by the computer model using the highest sea energy frequency component, c^max- 
This ensures that the linearly interpolated sea evolves in a smooth and continuous fashion. 

Its disadvantage is that it has poorer low frequency resolution of the sea wave energy than its 
SWS counterpart. This is demonstrated by a typical run case in Table 2.4 where the lowest 10 
frequencies defined are listed for both methods. As one can readily calculate from integrating 
Equation 2.7 with respect to a;, for a sea severity of 2.0 m significant wave height, 67.7% of the 
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Component 


Frequency, cjj 
SWS Method 


FFT Method 


1 


0.347 


0.347 


2 


0.381 


0.491 


3 


0.414 


0.601 


4 


0.448 


0.694 


5 


0.482 


0.776 


6 


0.515 


0.850 


7 


0.549 


0.918 


8 


0.583 


0.981 


9 


0.617 


1.041 


10 


0.650 


1.097 


Notes: 


1. FFT case: 
components 


Lj E [0.347, 22.2] with 4096 frequency 


2. SWS case: 
components 


u G [0.347,24.0] with 700 frequency 



Table 2.4: Comparison of Fi'equency Resolution for the SWS and FFT Methods 



sea wave energy lies below a frequency of 1.097 rad/sec. Thus the FFT method suffers in that 
it places most of the frequency components at high frequencies, where little of the wave energy 
resides. This undesirable result is a consequence of choosing a b space meshing with constant 
A6, which is necessary to compute FFT’s in the b and x domains. A constant Ab meshing 
produces a non-constant Aa; meshing with Acjj = \/27r^6„nn (\/i- s/i - l). 



2.4 Multidirectional Modeling Methods 

The current model improves upon the previous model by upgrading the two-dimensional unidi- 
rectional wave field discussed in Section 2.3 to a three-dimensional multidirectional wave field 
with the choice of any arbitrary tow angle in the wave field. 

The sea domain variables of interest for the multidirectional wave field are 

• sea elevation, h 

• horizontal fluid velocity (x-direction), 

• horizontal fluid velocity (y-direction), v 
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vertical fluid velocity, w. 



Two methods were considered in arriving at the sea domain properties of the multidirectional 
wave field, each with certain advantages and disadvantages. The first method will be called 
the Complete Wave Field (CWF) method, and the second will be called the Partial Wave Field 
(PWF) method. 

2.4.1 Complete Wave Field (CWF) Method 

The Complete Wave Field (CWF) method involves the construction of a three-dimensional wave 
field using a spectral energy model that is completely described in both the x and y directions 
in the spatial frequency domain. This is accomplished in a manner similar to that discussed 
in Section 2.3.2, except that the equations now will be a function of both x and y and will 
separate the horizontal velocity into x and y components. 

As previously mentioned in Equation 2.1, the sea elevation can be obtained by a superpo- 
sition of sine waves having a random phase and an amplitude prescribed by Equation 2.2. In 
like manner, the fluid velocities can be calculated by the superposition method. The formulas 
for /i, w, and w in an earth-fixed reference frame are provided in Equations 2.39 - 2.42. 



N N 



h{x', y', t) = cos{k^.{u;^.)x' + ky^ {u>y^)y' - ujjkt + <j)jk) (2-39) 



N N 




(2.40) 



X COs{kxj + ky^ {'■^yk)y' - >^jkt + 4>jk) 



N N 




(2.41) 



X COS(/Ca; . (WiJ x' + ky^ (cjyj I/' ~ LOjkt + <l>jk) 
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N N 



W 



Xj 5 ^Vk) siii(A:xj {^Xj ) 3 : + i^vk) y ^jki + ^jfc) (2.42) 

j=l k=l 



where 



Ojk = tan ^ 



Ayon . 



(2.43) 



^jk 



= \f^j 






(2.44) 



By invoking a coordinate transformation and the dispersion relation from Table 2.1, one 
arrives at the forms appropriate to the coordinate system moving with the antenna. These are 
presented in Equations 2.45 - 2.48. 



N N 

h{x,y,t) 

j=l k=l 



(2.45) 



X cos 



^Xj i^Xj ) 3 : + ky^, (Wyj. ) y (l^jk 



cos {0jk) 



+ ^jk 



N N 



U 



(x,y,t) = OJjkCOs(0jk) 

j=l k=l 



(2.46) 



X cos 



kx, (u’xj)x + ky^ (ojyJ y - {Ujk ~ 



cos {0jk) 



)t + ^jk 



N N 



’{x, y, t) = J2Y1 sin {0jk) 

j=i fc=i 



(2.47) 



X cos 



ujj^U cos(6>jfc). 



kxj{x>Xj) ^ ~k kyi^ {Uyi^) y i^jk ^ 4>jk 
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Bretschneider Energy Density Spectrum 




Figure 2-8: Complete Wave Field 3-D Energy Spectrum (quadrants 1 and 4 shown) 



w{x,y,t) 



N N 

^jki^Xj J ^yk ) ^jk 

j=\ k=l 



(2.48) 



X sm 



We again take advantage of the fact that Equations 2.45 - 2.48 are a Fourier series, and 
thus the functions /i,u,i;,and w are the inverse Fourier transform of some expression. Since 
the calculations to determine the sea domain properties occur at a fixed point in time, the 
transform variables are not uj and rather they become b and x. Because the two-dimensional 
FFT algorithm used for the model assumed the presence of frequency components in all four 
quadrants surrounding the origin in the spatial frequency domain, it was necessary to construct 
an energy density spectrum in three dimensions. An example of a three-dimensional energy 
density spectrum, in spatial frequency coordinates, which is used by the sea domain properties 
equations is shown in Figure 2-8. In order to retain the same amount of energy (i.e. area under 
the spectral energy density curve), it is necessary that the formula for the component amplitudes 
be adjusted from that given in Equation 2.2. In addition, it is convenient to incorporate the 
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random phase, into a complex component amplitude, cjk- With these changes in mind, the 
formulas providing the sea domain properties are given in Equations 2.49 - 2.54. 



h{x,y) 



N N 



j = l k = l 



(2.49) 



X exp 



i{27rbxjX -h + {27rbjkU cos 



(0jk) - \/2Tr gbjk) O] } 



N N 

u{x,y) = ^^Re|cjfc(6i^,6j,J\/2^^^cos(6>jfc) (2.50) 

j=i fc=i 

X exp |z(27r6i^a: + 2nby^y + {2nbjkU cos {6jk) - y/2'ngbjk) ^)j | 



N N 



’{x,y) = ^^Re{cjfc(5i,,5j,Jv'^7r^sjn(0jfc) (2.51) 

X exp ^i{2irbx.x + 2irby^y + {2irbjkU cos (djk) - >/2TTgbjk) O] | 



j=l k=l 



N N 

w{x,y) = ^^Re{-iCjfc(6a;^,6yJ \/27r^ (2.52) 

j=i fc=i 

X exp ^i{2TTb:cjX + 2nby^y + {2nbjkU cos {6jk) - \/2TTgbjk) O] | 



Rc{cjk{bxj,by^)} ^ ^bcos{(f>jk) 



(2.53) 



Im {cjfc(5xj,5yj} = ^^sin(0jfc) 



(2.54) 



where 



bjk = 



(2.55) 
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Then the Fourier transform pairs become 



H{b^,by) 




h{x,y) 


(2.56) 


s/2ngbjk cos {Qjk) H{bx, by) 




u{x,y) 


(2.57) 


^/2-Kgbjk sin {Qjk) H{bx, by) 




v{x,y) 


(2.58) 


-i^y2TTgbjkH{bx,by) 




w{x, y) 


(2.59) 



where 

N N 

H{b^, = (^jk{bxj,by^) exp (2TvbjkU cos {9jk) - y/^ngbjk^^ (2.60) 

j=l k=l 



H ^y) li{x, y) 



(2.61) 



is defined to mean 

N N 

H{bj,,by) = ^^/2.(x-j,yfc)exp(-z27r6^^x-i27rfe^^ (2.62) 

i=i k=i 
^ N N 

h{x,y) = — ^^//(6:i,^,fc2;Jexp(z27r6:i,^-x + z27r6y^ (2.63) 

^ j=i fc=i 

By setting the complex component amplitudes, Cjk, to obey the relations 



Cjk{ bxj^ by^) — Ci{bxj,byj^) (2.64) 

Cjk{-bxi,by^) = Ci(bx^,-byJ (2.65) 



that is, the frequency components of quadrants 2 and 3 are the complex conjugates of the 
frequency components of quadrants 4 and 1, respectively, then the imaginary parts of Equations 
2.49 - 2.52 are zero, and the Re{ } can be dropped from the equations. 

Finally, the wave field is constructed by performing an inverse FFT on the two-dimensional 
complex arrays, represented by the left side of Equations 2.56 - 2.59, which contain the Fourier 
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Three Dimensional Wave Field - CWF Method 



NOTE Only Portion of Wave Field Shovy/n 
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Figure 2-9: Wave Field Generated by CWF Method 

coefficients describing and w, A example of the wave field generated by the CWF 

method is shown in Figure 2-9. This computational step is not trivial and is the reason that 
tliis method was abandoned early in the project. The memory requirements and computation 
times necessary to complete the inverse FFT of these large arrays was not practical for desktop 
computer modeling. 

Wave field directionality has not been discussed up to this point. That topic has been 
reserved for Section 2.4.2, the Partial Wave Field method, which is the method that is imple- 
mented in the model. 

2.4.2 Partial Wave Field (PWF) Method 

The Partial Wave Field (PWF) method involves the construction of a three-dimensional wave 
field using a spectral energy model that is only partially described in both the x and y directions 
in the spatial frequency domain. In comparison to the CWF method, this partial description 
significantly reduces the amount of computation needed to build the modeled wave field. It 
reduces computation time by forming the multidirectional wave field from the proper combi- 
nation of several unidirectional waves as shown in Figure 2-10. In addition, the PWF method 
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Figure 2-10: PWF Unidirectional Wave Summation 

takes advantage of the fact that the only wave field information necessary to compute exposure 
statistics on the antenna is the information from the wave field that is directly in contact with 
the antenna. 

The first step of the PWF method is to select, from the available choices, the wave energy 
density spectrum to be modeled. The energy density spectrum dictates the behavior and 
magnitude of the modeled seas and its selection can significantly affect the exposure statistics. 
Guidance on proper selection of the energy density spectrum can be found in Section A. 2.1 and 
in numerous texts and papers. [4] [5] [7] [8] 

The next step of the PWF method is to construct the three dimensional energy density 
spectrum with NN discrete unidirectional spectra that are each propagating in a different 
direction. The model allows for the user to select any odd number of discrete unidirectional 
waves. The requirement for an odd number is driven by the assumption that one direction lies 
on the principal axis and the remaining directions are split evenly on either side of the principal 
axis. A larger number of directions will result in a more realistic three dimensional sea spectrum 
but at the expense of computation time. Referring again to Figure 2-10, one can see that as 
the number of unidirectional waves increases, the modeled seas more closely approach a realistic 
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Bretschneider Energy Density Spectmnn 




Figure 2-11: Partial Wave Field 3-D Energy Spectrum - No Directional Scaling (quadrants 1 
and 4 shown) 

three dimensional spectrum. Likewise, as the number of unidirectional waves decreases, the 
modeled seas move further away from a realistic three dimensional spectrum. The model 
constructs the three dimensional energy density spectrum by assigning a direction to each of 
the discrete unidirectional spectra, as shown in Figure 2-11. The example shown in Figure 
2-11 is for five propagation directions. For a given number of propagation directions, five in 
this example, the model assumes equal angles between each direction and then solves for the 
correct angle to assign to each direction. It should be noted that, for purposes of illustration, 
the magnitude of each of the unidirectional spectra is equal and that Figure 2-11 is similar 
to Figure 2-8 except that Figure 2-11 is only partially modeling the three dimensional energy 
density spectrum. 

It has already been pointed out the magnitude of each of the unidirectional spectra in Figure 
2-11 is equal. This implies that the wave field has no directional spreading characteristic; that 
is, energy is being propagated equally in all directions. In reality, wind-generated waves tend 
to propagate in the same general direction as the wind. Typically, the energy is spread over 
several directions, but the majority is in the direction of the wind. This spreading of energy 
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can be characterized by what is known cis a directional spreading function. The directional 
spreading function assigns a fraction of the total wave energy to each particular direction, with 
that fraction being determined by the shape of the spreading model selected. To generate 
a directional spectrum, one simply multiplies the energy density spectrum by the directional 
spreading function 



5 (a;, 9) = S {lv) D (0, a;) (2.66) 

where D (0,a;) represents the directional spreading as a function of both direction, 0, and fre- 
quency, u. Not unlike wave spectra, numerous directional spreading functions have been pro- 
posed by researchers attempting to characterize a particular sea wave environment. [4] [5] [6] [7] [8] 
Although the spreading functions implemented in the model are simple models that are a func- 
tion of direction only, spreading models exist that represent the spreading as a function of both 
direction and frequency. In these frequency dependent spreading functions, the low frequency 
components are assigned a very narrow spreading while the high frequency components are as- 
signed a wider spreading. [5] The spreading functions used in the model are the cosine square 
formula and the cosine fourth formula, given in Equations 2.67 and 2.68, respectively. 

D{e) = IcosHO), (2.67) 

7T Z Z 

D{e) = ( 2 . 68 ) 

A plot of the cosine square formula, shown in Figure 2-12, shows the behavior of the di- 
rectional spreading function as the propagation direction changes from 0 degrees, the principal 
axis of propagation, to ±90 degrees. As shown in Figure 2-12, the largest fraction of energy 
is along the principal axis while at ±90 degrees the fraction tapers to zero. To account for 
the total energy in the spreading spectrum, the actual fraction assigned to each unidirectional 
spectrum is the integral of the directional spreading function centered on each direction with 
integration limits being equal distances on either side of the propagation direction. Again, 
using the five propagation direction example, the propagation directions and their respective 
integration limits are shown in Figure 2-13. 

Applying the cosine square spreading function to the three dimensional energy density spec- 
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(^/^)cos^0 Directional Spreading Function 




Figure 2-12: Cosine Square Directional Spreading Function 



(^/ )cos^0 Directional Spreading Function 

71 




Figure 2-13: Cosine Square Directional Spreading Function for the Five Directions Example 
Problem 
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Brelschneider Energy Density Spectrum 




Figure 2-14: Partial Wave Field 3-D Energy Spectrum - With Directional Scaling (quadrants 1 
and 4 shown) 

trum in Figure 2-11, you obtain a three dimensional energy density spectrum with directional 
characteristics, shown in Figure 2-14. Now, comparing Figure 2-14 with Figure 2-11, you can 
see that the magnitudes of the unidirectional spectra in Figure 2-14 are not equal, but are scaled 
according to the directional spreading function. In addition, as prescribed by the spreading 
function, the principal axis of wave propagation, 0 degrees, contains the largest fraction of 
energy. 

The next step in the PWF method is to apply the FFT method, described in Section 2.3.2, 
to each of the scaled unidirectional wave spectra to solve for the spatial sea domain properties 
of each unidirectional wave. Now the individual waves must be properly summed to determine 
the three dimensional wave properties along the antenna. The sea domain properties of the 
PWF method are somewhat different than those found by the CWF method. In the CWF 
method, the goal was to solve for and w in the entire wave field. With the PWF 

method, the only information needed is the information along the antenna. To simplify the 
numerical modeling, the velocities of interest were found relative to the antenna. So rather 
than solving for the x and y velocities in the entire field, the model solves for the velocities 
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that act longitudinally and transversely on the antenna, as well as the vertical velocity. The 
summation is accomplished by Equations 2.69 - 2.72. 

NN 

{^) = X] fh (5) (2.69) 



NN 

^longitudinal (*^) ^ E Ui (x) cos (tp) (2.70) 

i=l 



NN 

^transverse (^) = E Ui{x) sin (tp) (2.71) 

i=l 



NN 

w{x) = Wi (x) 
i=l 



(2.72) 



where 



+ 0.5j Ascos(V^) (2.73) 

and 

NN is the number of propagation directions being summed 

hi (x) is the wave elevation of the individual unidirectional wave, at the location 5, in the 
unidirectional wave field that is acting on the antenna location .r, as shown in Figure 2-15 
tii (x) is the horizontal velocity of the individual unidirectional wave, at the location x, in 
the unidirectional wave field that is acting on the antenna location .r, as shown in Figure 2-15 
xjj is the angle between the individual unidirectional wave and the antenna, as shown in 
Figure 2-15 

lUi (x) is the vertical velocity of the individual unidirectional wave, at the location x, in the 
unidirectional wave field that is acting on the antenna location x, as shown in Figure 2-15 
Nsca is the number of spatial divisions in the unidirectional wave field 
As is the spatial increment on the antenna as well as the spatial frequency increment of the 
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Figure 2-15; Geometry of the Wave Summation 

unidirectional wave spectrum 

N is the number of spatial divisions on the antenna 
L is the modeled length of the antenna 

In Equations 2.69 - 2.73, it is assumed that the width of the unidirectional wave is infinite 
and that the surface properties along that infinite width do not vary. This assumption allows 
you to perform the summation shown in Figure 2-15. 

The origin of equation 2.73 is not obvious, but it represents the trigonometric and geometric 
relation between the individual wave and the antenna. Referring to Figure 2-15, the numbers 
represent the array storage locations of the particular variables with the wave properties num- 
bered i = 1 to Nsea and the antenna properties numbered j = 1 to W. The spacing between 
points, As, is the same for both the antenna and the wave with the length of the wave field 
being Lsea and the length of the antenna being L. In addition, it is assumed that the antenna 
and the wave cross at their midpoints. The position, x, can be found from 



(x, 



crossing point 



- X, 



origin 



) (^crossing point COS (^) 



(2.74) 
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where 



^crossing point 


= (number of whole and partial mesh spaces on wave field between 




crossing point and origin) x (As) 

= ('^“- 0 . 5 ) A. 


(2.75) 


^origin 


= 0.0 


(2.76) 


^crossing point 


= (number of whole and partial mesh spaces on antenna between 






crossing point and origin) x (As) 

= 0 . 5 ) As 


(2.77) 


X 


= (number of whole and partial mesh spaces on antenna between 






X and origin) x (As) 


(2.78) 




= U - 1) As 


(2.79) 



combining equations 2.74 - 2.79 you obtain 

X = ~ 0.5^ As — ~ j + 0.5^ As cos 

where j represents the array storage location nearest to x, given by 



(2.80) 



i = |(7V-l) + l 



(2.81) 



when Equation 2.81 is combined with Equation 2.80, you obtain 



X = 




|(iV^l) + l 




As cos [xl)) 



(2.82) 



which is identical to Equation 2.73. 

The PWF modeling framework consists of the following steps 

1. select the wave energy density spectrum to be modeled 

2. construct the three dimensional energy density spectrum with NN discrete unidirectional 
spectra that are each propagating in a different direction 
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3. select the directional spreading function 

4. apply the directional spreading function to the discrete components of the three dimen- 
sional energy density spectrum 

5. perform an inverse FFT on the one-dimensional complex arrays, represented by the left 
side of Equations 2.25 - 2.27, which contain the Fourier coefficients describing and 
w of each discrete unidirectional spectrum 

6. properly sum the sea domain properties of each discrete unidirectional wave along the 
length of the antenna to generate a three dimensional wave field on the antenna 

7. time advance the Fourier coefficients describing /i,n, and w using Equation 2.28 

8. repeat steps 5-7 for duration of simulation. 

2.5 Verification of Model Seas 

The random multidirectional seaway was constructed by the PWF method as described in 
Section 2.4.2. By running the modeling code and sampling the wave field for several parameters, 
one can build a wave record tliat can be analyzed to determine if the resultant wave field 
correctly models the desired wave field. The modeled seas verification focused on tlie following 
four areas 

• significant wave height, Hg 

• encounter frequency, cdg 

• spectral analysis of a spatial wave record 

• spectral analysis of a temporal wave record. 

2.5.1 Significant Wave Height 

By sampling the wave record for sea surface elevation, one can determine the significant wave 
height of the sample, as defined in Section 2.2. [7] This can then be compared to the significant 
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Sample 


Mean(/?i) 

(m) 


Spatial 

a 

(m) 


Hs 

(m) 


Temporal 

Mean (/i A a 

(m) (m) 


Hs 

(m) 


1 


-0.054G 


0.7177 


2.8779 


0.0102 


0.7784 


3.1215 


2 


-0.0120 


0.7455 


2.9893 


0.0032 


0.6950 


2.7869 


3 


0.0024 


0.8389 


3.3638 


-0.0009 


0.7893 


3.1652 


4 


O.OIGO 


0.8890 


3.5647 


0.0085 


0.7582 


3.0405 


5 


0.0287 


0.G33G 


2.5406 


0.0231 


0.7017 


2.8137 


G 


-0.0119 


0.8843 


3.5459 


-0.0045 


0.7159 


2.8707 


7 


O.OOGG 


0.7410 


2.9713 


0.0057 


0.7528 


3.0186 


8 


-0.0844 


0.7G3G 


3.0621 


-0.0171 


0.7140 


2.8632 


9 


-0.03G2 


0.5490 


2.2016 


-0.0185 


0.7412 


2.9723 


10 


-0.0391 


0.85G9 


3.4360 


0.0092 


0.7373 


2.9567 


11 


0.0G32 


0.7047 


2.8258 


-0.0182 


0.7215 


2.8934 


12 


0.0050 


0.6758 


2.7100 


-0.0061 


0.7636 


3.0619 


Averages: 


-0.0097 


0.7500 


3.0075 


-0.0005 


0.7391 


2.9637 



Table 2.5: Simulated Significant Wave Heights for a Specified Sea Severity of 3.0 m 



wave height which was specified in the input file. The measured significant wave height can be 
found by Equations 2.83 and 2.84, 



a 



n 



z=l 

Hs = 4.01 



(2.83) 

(2.84) 



wliere 

Xi is the sampled wave amplitude, measured from the mean free-surface 
a is the standard deviation of the sampled wave amplitudes 
II s is the measured significant wave height of the wave record. 

(The “ ” symbol denotes that the quantity is computed from a sampled wave record, and 

therefore may deviate somewhat from the theoretical value.) 

A sea severity of 3.0 m significant wave height was simulated by the PWF method. Twelve 
wave amplitude records were analyzed by sampling the seas spatially at 0.0625 m intervals 
and temporally at 0.021 sec intervals. The results are summarized in Table 2.5. An excellent 
agreement to the specified significant wave height is achieved with an average deviation in Hg 
of 0.25% and 1.20% for the spatial and temporal analysis, respectively. 
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Encounter Frequency for Varying Tow Speeds 




Figure 2-16: Encounter Frequency for Varying Tow Speeds 



2.5.2 Encounter Frequency 

Determination of the encounter frequency presents a formidable task for a random multidirec- 
tional seaway. To simplify the task, the analysis of encounter frequency was done using the 
single sine wave spectrum as the input spectrum. The method used to verify the encounter 
frequency invoked the superposition principle: the sea waves generated by the PWF method are 
ultimately a superposition of individual sine waves, thus the verification of encounter frequency 
was restricted to seas represented by a single sine wave. The goal of this analysis is to verify 
that the antenna sees the correct “encounter frequency,” i.e. an effective change of the wave 
frequency as the antenna is towed in a seaway. The encounter frequency, cjg, is given by 



— ^0 



cJqU cos {9) 
9 



(2.85) 



where 

loq is the wave frequency in a stationary coordinate system 
U is the antenna tow speed 

9 is the tow angle with respect to the principal axis of wave propagation 



45 



Spatial Power Spectral Density Analysis 



1 Propagation Direction 




Figure 2-17: Spectral Analysis of Spatial Wave Data - One Propagation Direction 
g is the acceleration due to gravity. 

Thus for towing into head seas (cos (180®) is negative), the encounter frequency goes up, 
and for towing with following seas (cos(0®) is positive), the encounter frequency goes down. 
The encounter frequency was measured for a single sine wave of frequency cjq = TO rad/ sec at 
tow speeds from 0.0 mj sec to 5.0 m/ sec at a tow angle of 45®. The resulting plot of vs. U 
is shown in Figure 2-16. The plot of uJe vs. U is expected to be linear with a slope of 
(see Equation 2.85), which for cjq = 1.0 rad/ sec and 6 = 45® has a value of -0.0721 m“^. The 
measured slope was -0.0767 which shows the proper dependence of encounter frequency 
on speed has been captured by the model. 

2.5.3 Spatial Spectral Analysis 

By sampling the wave record spatially for sea surface elevation, one can determine the spatial 
power spectral density curve of the modeled wave field. By comparing the resultant spectral 
density curve to the desired spectral density curve, a verification of the model’s wave field 
generation method can be conducted. Figure 2-17 shows the spectral analysis of spatial wave 
data for a wave field constructed with only one propagation direction. The resultant energy 
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spatial Power Spectral Density Analysis 
3 Propagation Directions 




Figure 2-18: Spectral Analysis of Spatial Wave Data - Three Propagation Directions 

density spectrum closely resembles the desired input spectrum and there is good agreement 
amongst the values representing the area under the input spectrum, the area under the output 
spectrum, and the mean square wave elevation of the modeled wave field. 

An interesting phenomena occurs when the number of propagation directions is increased 
to three. As shown in Figure 2-18, there is an apparent increase in energy density at the lower 
spatial frequencies while the opposite occurs at the higher spatial frequencies. There is still 
good agreement amongst the values representing the area under the input spectrum, the area 
under the output spectrum, and the mean square wave elevation of the modeled wave field, 
but there is an obvious shift in energy to lower frequencies. This shift in energy to lower 
frequencies is attributed to the wave summation method listed in Equation 2.73 and shown in 
Figure 2-15. As the individual wave is added, it is being ’’stretched” in length due to the angle 
between the antenna and the wave itself. For the extreme case where this angle is 90°, only one 
value of wave elevation from the individual wave would be applied to the entire length of the 
antenna. In this case, a single point value is stretched along the entire length of the antenna. 
In the opposite extreme case where the angle is 0°, there would be a point to point match in 
the individual wave and the antenna with no stretching, as shown in Figure 2-17. The same 
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Spatial Power Spectral Density Analysis 
5 Propagation Directions 




Figure 2-19: Spectral Analysis of Spatial Wave Data - Five Propagation Directions 

concept applies to angles in between the extreme cases as a small portion of the wave field is 
stretched to fit the antenna. 

The stretching causes a decrease in spatial frequency as the apparent wavelength of the 
individual wave is increased, as shown in Figure 2-20. Since the individual spatial frequency 
decreases, the total energy of the wave field along the antenna shifts to lower frequencies. It 
should also be pointed out that the shifts in energy magnitude occur at values that correspond 
to the cosine of the angle between the individual waves and the antenna. In the case of three 
propagation directions, shown in Figure 2-18, the waves are propagating at -60°, 0°, and -hG0°. 
The shift in energy magnitude occurs at 0.5 which is the same value as the width of the 
input spectrum multiplied by cos (60°) . In the case of five propagation directions, shown in 
Figure 2-19, the waves are propagating at -72°, -36°, 0°, +36°, and +72°. The shifts in energy 
magnitude occur at 0.3 and 0.8 which are the same values as the width of the input 
spectrum multiplied by cos (36°) and cos (72°), respectively. 
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Figure 2-20: Spatial Frequency Stretching 
2.5.4 Temporal Spectral Analysis 

By sampling the wave record temporally for sea surface elevation, one can determine the tem- 
poral power spectral density curve of the modeled wave field. By comparing the resultant 
spectral density curve to the desired spectral density curve, a verification of the model’s wave 
field generation method can be conducted. Figure 2-21 shows the temporal spectral analysis 
of a model generated wave field. The resultant energy density spectrum closely resembles the 
desired input spectrum and there is good agreement amongst the values representing the area 
under the input spectrum, the area under the output spectrum, and the mean square wave 
elevation of the modeled wave field. 
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Temporal Power Spectral Density Analysis 




Figure 2-21: Spectral Analysis of Temporal Wave Data 
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Chapter 3 



Antenna Model Development 



3.1 Coordinate System Description 

The coordinate system adopted for the model is a three-dimensional system whose origin trans- 
lates in the antenna tow direction at the antenna tow speed [/, which is defined to be positive. 
The origin of the antenna is defined to be the trailing point of the antenna. The waves travel 
in the positive x direction, all at varying angles, but always from left to right in the positive x 
direction. Figure 3-1 depicts the coordinate system in the x-z plane as well as the basic antenna 
configuration. The antenna tow angle, 0, is the angle measured from the antenna to the prin- 
cipal axis of wave propagation. It is defined to be 0^ when pointing in the positive x direction 
and 180° when pointing in the negative x direction. The model allows the user to select any 
tow angle between and including 0° and 180°. Due to port and starboard symmetry of the 
anntena, the behavior of the antenna for the angles 180° to 360° is identical to the behavior 
of the antenna for the angles 0° to 180°. Figure 3-2 depicts the coordinate system in the x-y 
plane as well as the antenna tow direction. 

The equations listed in Table 2.1 apply to a stationary coordinate system. The transforma- 
tion to a moving coordinate system employs the following relationship 



— X -{-Ut 



(3.1) 



where 
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Direction of 




Figure 3-1: Coordinate System: X-Z Plane 




Following 

Seas 



Figure 3-2: Coordinate System: X-Y Plane 
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x' is the earth-fixed, stationary, coordinate system 

X is the moving coordinate system whose origin translates with velocity U , as shown in 
Figures 3-1 and 3-2. 

3.2 Governing Vertical Equation of Motion 

The vertical forces acting on the antenna are given in the right hand side of Equation 3.2 



where 

m is the mass per unit length of the antenna, assumed uniform 
V is antenna vertical velocity, which is equal to ^ 

T{x^t) is the tension along the length of the antenna, arising from skin frictional drag 
between the submerged portion of the antenna and the water 

q is the z coordinate of the antenna central axis, as shown in Figure 3-3 
Cd is the normal drag coefficient of a cylindrical antenna section as it moves in the z 
direction through the water 

w is the fluid vertical velocity 

b{x^ i) is the net buoyancy force acting on the antenna, which is the difference between the 
downward force due to its weight and the upward force due to the displaced water. 

Cm A is the coefficient of added mass for a cylindrical section in a flow field normal to its 
long axis 

Pwater i^ density of water (for seawater of 15° C, py^ater ~ 1025.9 kg/m^) 
r is the antenna radius 

El is the product of the antenna cross sectional modulus of elasticity and area moment of 
inertia, thus a measure of lateral stiffness. 

The first term on the right hand side of Equation 3.2 represents a restoring force charac- 
teristic of a string displaced from its equilibrium condition. The tension in the antenna is 
calculated from the time varying wetted surface area along the length of the antenna. A tan- 
gential drag coefficient of Ct = 0.0035 was used based upon Reynold’s numbers which ranged 
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Figure 3-3: Antenna and Sea Surface Elevation Variables 

from 32,940 to 592,900 for speed /diameter coiiibiiiatioiis of 3 knots/1 inch to 9 knots/6 inches, 
respectively. [1] The second term represents the normal drag force a cylindrical element would 
experience as it is moved through the fluid with its long axis perpendicular to the direction of 
motion. A value of Co = 1.0 was used, which assumes an infinite fluid domain. The fourth 
term represents the added mass effect experienced by a cylindrical object as the fluid around it 
is accelerating normal to the long axis of the cylinder. A value of Cm a = 1.0 was used, which 
is consistent with standard fluid dynamic texts assuming an infinite fluid domain. [6] The final 
term represents the restoring force due to the stiffness of the material, and its distribution with 
respect to the neutral axis. Initial design options have the lateral stiffness ranging from 300 to 
1400 Pa-w?. 

It is recognized that some of these forces have been approximated, namely in that the values 
quoted above are for Co and Cm A hi an infinite fluid domain. The antenna operates in close 
proximity to the free surface, and occasionally portions rise completely above the surface of 
the water. A rigorous treatment of the behavior of Co and Cm A near the free surface was not 
attempted. Rather, a straightforward linear approximation was employed which varied the 
coefficients Co and Cm A from their infinite fluid domain values when the antenna axis was at 
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Conductor: Antenna Element: 

Isometric View Section View 




Figure 3-4: Antenna Element Description 

a depth of r, just submerged, to a value of zero when the antenna axis was at a height of r, 
just breaking out of the water. This dependence of Co and Cm A on relative antenna depth, y, 
is shown in Equations 3.3 and 3.4, which only apply to the partially submerged case. 



CD{y) = + ^) 


(3.3) 


y 1 

CMA{y) = C'A/-4/n/ini(e(-^ + 2 ) 


(3.4) 


3.3 Antenna Modeling Methods 




3.3.1 Discretization 





The antenna, of length L, is subdivided first into lengths of the physical radio antenna elements. 
These antenna elements, designed at Lincoln Laboratory, are called “jelly roll” antennas because 
they are a conducting sheet rolled up about 1^ times. A graphical description of the antenna 
element is shown in Figure 3-4. The element length is designated Leiem- This first level of 
subdivision is not only logical, but necessary for the exposure statistics calculated by the model. 
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Antenna 

Length 

L 



Element 
Length ^ 
L. 



■^elem 



^elem 



Element Length 



■^elem 



^ ^elcm ^elem 



Lelem = M As 



Mesh 
Length ^ 
As 



M = 6 



Figure 3-5: Antenna and Antenna Element Description 

As it turns out, this type of antenna element performs well as long as it has some small amount, 
on the order of 1 mm, of its outer housing diameter exposed along its entire length. But if 
the element becomes covered by even a small amount of water at any point along its length, 
antenna performance drops dramatically. Thus it is necessary to monitor numerous points 
within each element to determine whether a given element is exposed or not. This leads to 
the final level of antenna subdivision: each element is meshed into some integer number M of 
sections with length As. The subdivision process is shown in Figure 3-5. 

An interesting feature of the models which calculate the sea elevation and fluid velocities 
by Fast Fourier Transforms (FFTs) is that the frequency meshing (Aa;) of the spectral energy 
density function S{u)^ is now related to the spatial meshing (Ax) of the sea, and thus the 
spatial meshing of the antenna (As). This means that there is an extra constraint to consider 
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in selecting the finest antenna mesh size As; or conversely, given a mesh size As, constraints 
are imposed on the Au; meshing of the spectral energy density function Some freedom is 

regained if one decouples the length of the antenna and the length of the sea simulated. This 
decoupling is implemented in the model, so that for a typical run the antenna may be 25 m 
long with a mesh size of 0.0625 m, while the length of sea modeled is 512 m long with a mesh 
size of 0.0625 m. The key point is that the mesh length between the antenna (As) and the sea 
(Ax) must match. 

3.3.2 Integration Method 

The time marching integration method implemented is an explicit forward Euler scheme. The 
following discussion of of the explicit forward Euler scheme comes from Ulrich. [9] Such schemes 
are characterized by using information at time t to calculate functions at time ^ + A^. A simple 
example of how the forward Euler method is used is given below. Consider the following sample 
numerical problem. 

Solve 



du 

dt 



u{0) 



1 



The time derivative is approximated using the Euler method as 



du 

dt 



t=t^ 



jjn+l _ ijn 



which for constant time step A^ becomes 



du 

dt 






_ ijn 

A^ 



having adopted the notation of Celia and Gray [2], where 



(3.5) 

(3.6) 



(3.7) 



(3.8) 



f/" = u(r) 



(3.9a) 
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and is the time step. Then the solution proceeds step by step, beginning with the initial 
value of u given in Equation 3.6 



du 

dt 



= u(0)2 = 1 



t=o 



(3.10) 



then solving Equation 3.8 for we obtain 



du 



=U^ + At — 
dt 



(3.11) 



t=o 



With the value of u now known at time t (C/^), we can calculate 



du 

dt 



^2 



t=i 



= {U^) 



(3.12) 



and 



(3.13) 

t=0 

In like manner, we move forward in time to calculate any desired f/^. Numerical experi- 
mentation was used to determine the appropriate time step size which produced acceptable 
results. 

The spatial derivatives used in the model are classified as centered finite-divided- 

difference formulas. [3] Thus the antenna algorithm can be classified as a forward time, centered 
space (ETCS) method. 

3.3.3 Boundary Conditions 

The boundary conditions imposed on the antenna endpoints are given in Equations 3.14 and 
3.15 as 



U^ = U^+At ^ 
dt 



9(0) 


= h{0) 


(3.14) 


q{l) 


= h{L) 


(3.15) 
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That is, the antenna follows the sea surface at its endpoints. This is envisioned to be accom- 
plished at the trailing end by a drogue having positive buoyancy, and at the tow point by one 
of several options. Tliese options include 



• a dynamic lifting body 

• a surface following buoy 

• a leader of positive buoyancy 

• some combination of the above features. 

3.3.4 Initial Conditions 

The antenna simulation begins at i = 0 with 



q{x) = h{x) 

dq{x) dh{x) 

dt dt 



(3.16) 



(3.17) 



that is, the antenna is placed at the sea surface, and given an initial vertical velocity matching 
that of the sea surface. 
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Chapter 4 



Modeling Results 



4.1 Exposure Statistics 

4.1.1 Definition of Exposure 

The primary goal of the modeling effort was to predict exposure statistics as various design 
and operational parameters were varied; therefore, it is necessary to define exposure before 
proceeding on to the exposure statistics results. The overall exposure of the antenna is driven 
by the exposure of individual antenna elements, so the determination of an exposed, or receiving, 
antenna element must be addressed. The experimentally observed performance characteristics 
of the “jelly roll” antenna elements led to the following definition of exposure. It was found that 
for a level element in calm water, the antenna gain is stable right up to the point at which the 
outer skin becomes completely submerged. It was also found that if any portion of the outer 
skin was completely covered by water then the gain of the element suffered drastically. Thus an 
exposed element was defined to be one which was exposed by some small threshold height, set 
to four percent of Diameter for all simulations, at every point along the length of the element. 
In a discrete computer model some approximation is necessary with regard to every point along 
the length of the element. Thus each antenna element was divided into smaller segments for 
analysis. Typical designs considered had an element length of 0.5 m and were further divided, 
or oversampled, at M = 8. 
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4.1.2 Statistics Calculated by the Model 



The computer model was programmed to track three key statistics: 

• average exposure - the average fraction of antenna elements exposed throughout the entire 
duration of the simulation 

• threshold met - the average fraction of the time that some specified number of antenna 
elements, or more, were exposed 

• average counterdetection length - the average length of the antenna where there is no 
water touching the antenna. 

In addition, a record is made at each sample time of the condition of every antenna element, 
such that when the simulation is finished, an exposure matrix is available for further post- 
processing. Examples of information which could be extracted from the exposure matrix 
include 

• average drop-out time of a specific element 

• average spatial length of wave wash-overs 

• histograms and standard deviations of the number of exposed elements 

• exposure statistics for an array not fully populated (i.e., one could examine only the 
columns which had antenna elements installed) . 

Five output files are generated for each run of the computer model. These output files are 
described in detail in Appendix A. The input parameters which are specified for each run of 
the computer model are also described in detail in Appendix A. 

In the exposure statistics results given below, each series of simulations is presented in the 
following manner: 

• a listing of the antenna and operational parameters input table used to generate the series 

• a plot of the average fraction of elements exposed 
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• a brief discussion of the exposure plot. 

It is important to realize that these plots represent averages in time and space. This is 
necessary in order to display the results as functions of key antenna design and operational 
parameters, which is most useful from a designer’s viewpoint. FVom the output files, it is also 
possible to analyze how the statistics vary as time progresses, or how the statistics vary along 
the array from the tow point to the trailing end. Although there are a large number of possible 
statistics that could be analyzed for an even larger number of design and operational param- 
eter combinations, only one statistic is presented here for seven selected design or operational 
parameters. 
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Parameter Description 



Duration 

^ el em 

Lant 

Lelem 

Diameter 

El 

Pant 

M 

U 

Tow Angle 
Hs 

Development 
Spectral Model 
Spreading Function 

^ freq 

NN 



60 sec 
73 

36.5 m 
0.5 m 

3.0 inches 
100 N-m^ 

408 kg/m^ 

8 

5.0 kts 

0 ^ 45 ^ 90 ^ 135 ^ 180 ° 

3.0 m 

1.0 

Bretschneider 

2 

COS 

213 

1,5,9,13,17,21 



Table 4.1: Series 1 Input Parameters 



4.2 Results 

4.2.1 Series 1; Statistics vs. Number of Propagation Directions 
Input 

The input parameters for Series 1: Statistics versus Number of Propagation Directions, are 

given in Table 4.1. This is a measure of the fineness of the angular sea description. The 
number of propagation directions used were 1, 5, 9, 13, 17, and 21. The tow angles used were 
0, 45, 90, 135, and 180 degrees. 
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Exposure vs Number of Propagation Directions 




Figure 4-1: Average Exposure vs. Number of Propagation Directions 



Results 

Figure 4-1 shows a very weak dependence of the average fraction of elements exposed on the 
number of propagation directions used by the model to simulate the multidirectional wave field. 
However, the plot shows a strong dependence of average exposure on tow angle. The insensi- 
tivity to the number of propagation directions in the sea description not only demonstrates the 
validity of the approach used, but also shows that computing time can be minimized by using 
just a modest number (3 to 9) of directions in the sea description. 
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Parameter 


Description 


Duration 


60 sec 


^elem 


73 


Lant 


36.5 m 


^elem 


0.5 m 


Diameter 


3.0 inches 


El 


100 N-m? 


Pant 


408 kg/m^ 


M 


8 


U 


5.0 kts 


Tow Angle 


0^30^60^90^ 1 20^ 150^ 1 80° 


Hs 


3.0 m 


Development 


1.0 


Spectral Model 


Bretschneider 


Spreading Function cos^ 


^ freq 


213 


NN 


3 



Table 4.2: Series 2 Input Parameters 



4.2.2 Series 2: Statistics vs. Tow Angle 
Input 

The input parameters for Series 2: Statistics versus Tow Angle, are given in Table 4.2. The 
tow angles used were 0, 30, 60, 90, 120, 150, and 180 degrees. 
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Exposure vs Tow Angle 




Figure 4-2: Average Exposure vs. Tow Angle 



Results 

Figure 4-2 shows a rather weak dependence of the average fraction of elements exposed on the 
tow angle from 0^-90° (following seas); however, there is a very strong dependence on tow angle 
from OO'^-ISO® (head seas). 
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Parameter 


Description 


Duration 


60 sec 


^elevn 


73 


L'ant 


36.5 m 


L'elem 


0.5 m 


Diameter 


3.0 inches 


El 


100 N-m^ 


Pant 


408 kg/m^ 


M 


8 


U 


5.0 kts 


Tow Angle 


0^45^90^135^180^ 


Hs 


0.5,1.0,1.5,2.0,2.5,3.0 m 


Development 


1.0 


Spectral Model 


Bretschneider 


Spreading Function 


2 

cos 


^ freq 


213 


NN 


3 



Table 4.3: Series 3 Input Parameters 



4.2.3 Series 3: Statistics vs. Sea Severity 
Input 

The input parameters for Series 3: Statistics versus Sea Severity, are given in Table 4.3. The 
significant wave heights used were 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 meters. The tow angles used 
were 0, 45, 90, 135, and 180 degrees. 
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Exposure vs Significant Wave Height 




Figure 4-3: Average Exposure vs. Sea Severity 



Results 



Figure 4-3 shows a rather weak dependence of the average fraction of elements exposed on the 
significant wave height at all tow angles. 
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Parameter 


Description 


Duration 


60 sec 


^elem 


73 


L/ant 


36.5 m 


Lfelem 


0.5 m 


Diameter 


3.0 inches 


El 


200 N^m? 


Pant 


408 kg/rn^ 


M 


8 


U 


3,5, 7,9 kts 


Tow Angle 


0°,45°,90°,135°,180° 


Hs 


3.0 m 


Development 


1.0 


Spectral Model 


Bretschneider 


Spj'eading Function cos^ 


^ freq 


213 


NN 


3 



Table 4.4: Series 4 Input Parameters 



4.2,4 Series 4: Statistics vs. Tow Speed 
Input 

The input parameters for Series 4: Statistics versus Tow Speed, are given in Table 4.4. The 
tow speeds used were 3, 5, 7, and 9 knots. The tow angles used were 0, 45, 90, 135, and 180 
degrees. 
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Exposure vs Tow Speed 




Figure 4-4: Average Exposure vs. Tow Speed 



Results 

Figure 4-4 demonstrates the efTect that tow angle has on other parameters. In this instance, 
increasing tow speed for tow angles greater than 90^ (head seas) causes a decrease in exposure 
fraction. While, for tow angles less than 90° (following seas), increasing tow speed causes an 
increase in exposure fraction. 
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Parameter 


Description 


Duration 


GO sec 


^elem 


10,30,50,120,200 


I^ant 


5,15,25,60,100 m 


^elem 


0.5 m 


Diameter 


3.0 inches 


El 


200 N-rri^ 


Pant 


408 kg/m^ 


M 


8 


U 


5.0 kts 


Tow Angle 


0^45^90^135^180^ 


Hs 


3.0 m 


Development 


1.0 


Spectral Model 


Bretschneider 


Spreading Function cos^ 


Efreq 


213 


NN 


3 



Table 4.5: Series 5 Input Parameters 



4.2.5 Series 5: Statistics vs. Antenna Length 
Input 

The input parameters for Series 5: Statistics versus Antenna Length, are given in Table 4.5. 
The antenna lengths used were 5, 15, 25, 60, and 100 meters. The tow angles used were 0, 45, 
90, 135, and 180 degrees. 



71 



1 



Exposure vs Antenna Length 




Figure 4-5: Average Exposure vs. Antenna Length 



Results 

Figure 4-5 shows a very strong dependence of the average fraction of elements exposed on the 
antenna length for lengths less than 25 meters. For lengths greater than 25 meters, there is 
little or no dependence on antenna length. 
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Parameter 


Description 


Duration 


60 sec 


^elem 


73 


Lant 


36.5 m 


Lelem 


0.5 m 


Diameter 


3.0,6. 0 inches 


El 


200 N-m'^ 


Pant 


408 kg/m^ 


M 


8 


U 


5.0 kts 


Tow Angle 


0°, 45°, 90°, 135°, 180° 


Hs 


3.0 m 


Development 


1.0 


Spectral Model 


Bretschneider 


Spreading Function cos^ 


^ freq 


213 


NN 


3 



Table 4.6: Series 6 Input Parameters 



4.2.6 Series 6: Statistics vs. Antenna Diameter 
Input 

The input parameters for Series 6: Statistics versus Antenna Diameter, are given in Table 4.6. 
The antenna diameters used were 3 and 6 inches. The tow angles used were 0, 45, 90, 135, 
and 180 degrees. 
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Exposure vs Antenna Diameter 




Figure 4-G: Average Exposure vs. Antenna Diameter 



Results 

Figure 4-6 also demonstrates the influence that tow angle has on amplifying the effects of design 
parameter changes. In all cases shown, the larger diameter antenna has an increased exposure 
fraction, but the change is most significant for head seas cases. 
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Parameter 


Description 


Duration 


60 sec 


^elem 


73 


I^ant 


36.5 m 


Lclem 


0.5 m 


Diameter 


3.0 inches 


El 


200 N-ni^ 


Pant 


200,300,400,500,600 kg/n? 


M 


8 


U 


5.0 kts 


Tow Angle 


0°,45°,90°,135°,180° 


Hs 


3.0 m 


Development 


1.0 


Spectral Model 


Bretschneider 


Spreading Function 


cos^ 


^ freq 


213 


NN 


3 



Table 4.7: Series 7 Input Parameters 



4.2.7 Series 7: Statistics vs. Antenna Density 
Input 

The input parameters for Series 7: Statistics versus Antenna Density, are given in Table 4.7. 
The antenna densities used were 200, 300, 400, 500, and 600 kg/m^. The tow angles used were 
0, 45, 90, 135, and 180 degrees. 
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Exposure vs Antenna Density 




Figure 4-7: Average Exposure vs. Antenna Density 



Results 

Figure 4-7 again demonstrates the amplifying effect that tow angle has on design parameter 
changes. In all cases shown, as you increase the antenna density there is a resulting decrease in 
exposure fraction; however, this decrease is amplified for tow angles in the head seas direction. 
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Chapter 5 



Future Efforts 



The combined Navy and DARPA Multi-element Buoyant Cable Antenna effort was funded 
in FY 00 for 3 years. The demonstration phase culminated with a successful testing effort in 
Hawaii in the Summer of 1999. Within FY 00, Lincoln Laboratory is scheduled to complete the 
development and design of an antenna which can stand up to the rigors of an at sea deployment. 
FY 01 will see the initial production of the deployable version with the goal of at sea submarine 
testing during FY 02. 

Work at MIT will continue to support the buoyant cable antenna effort. Another graduate 
student, under the supervision of Professor Jerome Milgrarn, will continue and extend the 
scope of the simulation effort. Areas which will be investigated to improve the capability of 
the computer model include 

• more detailed modeling of the hydrodynamics at the free surface 

• higher order integration schemes 

• addition of digital signal processors to improve computation speed 

• validation against tow tank data 

• validation against at-sea test data 

• effects of roll for revised antenna elements whose performance depends on roll angle. 
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Appendix A 



User’s Guide to Proteus 

A-1 Introduction 

A. 1,1 Contents of the User’s Guide 

The User’s Guide is broken down into three sections: 

• Introduction 

• Input File 

• Output Files, 

Section A.l describes the contents of the user’s guide and contains the procedures necessary 
to obtain the results discussed in Section A. 3. Section A. 2 describes the Proteus input file line 
by line and includes program features and limitations associated with each parameter. Section 
A. 3 briefly describes the contents of each of the five output files and the recommended methods 
of viewing or extracting the information. 

A. 1.2 Analysis Procedure 
Obtaining the Snapshot Output 

1. Ensure that the input file, proteus_inputdxt, is located in the same directory as the 
FORTRAN executable file, proteus.exe. 
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2. Modify the input file to reflect the conditions that the user wishes to analyze, using Section 
A. 2 as guidance for features and limitations. Ensure that the rnatlab_ output variable is 
set to .TRUE. 



3. Execute proteus.exe, 

4. The snapshot output file /' will be located in the same directory as the 

executable file. 



5. 



6 . 



Ensure that the output file (nin_Z(i//’_mat’) will be located in the same directory as 
the Matlab^ function file, proteus_ mat 
(R) 

At the Matlab^ command line type: clear; proteus_mat(Vun_id/ / _maV) , 



Obtaining the Parameters and Results Summary 

1. Ensure that the input file, proteus_input.txt^ is located in the same directory as the 
FORTRAN executable file, proteus.exe, 

2. Modify the input file to reflect the conditions that the user wishes to analyze, using Section 
A.2 as guidance for features and limitations. Ensure that the 7xtn_ info variable is set to 
.TRUE. 

3. Execute proteus.exe. 

4. The parameters and results summary output file {run_id / /^ _\nio.txV) will be located in 
the same directory as the executable file. 

5. The output file can be viewed with any standard text editor. 



Obtaining the Exposure Data Matrix 

1. Ensure that the input file, proteus_input.txf is located in the same directory as the 
FORTRAN executable file, proteus.exe. 

2. Modify the input file to reflect the conditions that the user wishes to analyze, using Section 
A.2 as guidance for features and limitations. Ensure that the exp_ data variable is set to 
.TRUE. 
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3. Execute proteus.exe. 



4. The exposure data matrix output file (ran_i^i//’_edm’) will be located in the same 
directory as the executable file. 



5. 



(r) 

The output file can be viewed with a user supplied Matlab^ 
application. 



file or with a spreadsheet 



Obtaining Wave Data Temporal Spectrum Analysis 

1. Ensure that the input file, proteus_input,txt^ is located in the same directory as the 
FORTRAN executable file, proteus.exe. 

2. Modify the input file to reflect the conditions that the user wishes to analyze, using Section 
A. 2 iis guidance for features and limitations. Ensure that the wave_temporal variable is 
set to .TRUE. 



3. Execute proteus.exe. 

4. The temporal wave data output file (run_ irf//’_mat2’) will be located in the same di- 
rectory as the executable file. 



5. Ensure that the output file (n/n_irf//’_mat2’) will be located in the same directory as 
the executable file, proteus_mat2.exe. 



6. Execute proteus_mat2.exe. When the program prompts the user for filename, enter 
7mn_ id/ /_ mat2. 



7. 



Ensure that the output file (ru7i_zrf//’_mat2_proc’) will be located in the same directory 
as the Matlab® function file, proteus_ mat2. 



8. At the Matlab 






command line type: 



clear; proteus_ mat2(^run_ id/ /_ mat2_proc^). 



Obtaining Wave Data Spatial Spectrum Analysis 

1. Ensure that the input file, protcus_input.txt., is located in the same directory as the 
FORTRAN executable file, proteus.exe. 
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2. Modify the input file to reflect the conditions that the user wishes to analyze, using Section 
A. 2 as guidance for features and limitations. Ensure that the wave_spatial variable is 
set to .TRUE. 



3. Execute proteus.exe. 

4. The spatial wave data output file {run_ zd//’_mat3’) will be located in the same directory 
as the executable file. 



5. Ensure that the output file (ran_zd//’_mat3’) will be located in the same directory as 
the executable file, proteus_mat3.exe. 



6. Execute proteus_mat3.exe. When the program prompts the user for filename, enter 
run_ id/ /_ mat3. 



7. 



Ensure that the output file (ran_Z(i//’_mat3_proc’) will be located in the same directory 
as the Matlab^ function file, proteus_mat3. 



8. 



At the Matlab 






command line type: 



clear; proteus_ mat3(’run_ id/ / _ mat3_ proc ’). 



A. 2 Proteus Input File 

A. 2.1 Description of Features and Limitations 

Header 

The input file for Proteus is shown in Section A. 2. 2. The first 7 lines are for informational 
purposes and have no effect on the program. However, note that all variables are given in SI 
units and that there must be 7 lines prior to the first variable defined. 



Run Identifier (run^id) 

The run_ id is an 8 character variable used to name the 5 output files. The 5 output files are 
named 

• run_id/ _ma,V 
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• n/72_ui//’_info.txt’ 

• run_ id/ /’_edm’ 

• nm_ u///’_mat2’ 

• 7'7m_zd//’_inat3\ 

(The / / represents FORTRAN90 string concatenation.) Output files are located in the local 

directory from which the program is run. A sample pattern for run_id is nm_ zd=T5May_0r. 

(r) 

Characters which Matlab^ would interpret as operators (/, +, *, etc.) should be avoided in 

run_id if one of the Matlab® output options is selected. 

FFT Power of 2 {nnf) 

The FFT algorithm is designed to o])erate on arrays which have an integer power of 2 elements. 
nnf sets this power of 2, making 






(A.l) 



where 

N _se a is the number of sea frequency components and spatial locations. 

The length of the sea surface modeled in the spatial domain is generally much larger than 
the antenna length requires; however, the speed advantage of the FFT method more than 
outweighs the penalty of carrying along the extra sea information. 

Number of Antenna Elements {num^of^elern) 

This is the number of physical antenna receiver elements in the array, assuming that 

• the array is fully populated 

• the elements have zero space between them. 

Antenna Element Length {element _ length) 

This is the length of each antenna element. 
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Antenna Oversample Factor (M) 

Each antenna element is broken down into an integer number, M, of sections. These sections 
correspond to the antenna spatial mesh size As. By enforcing the condition that As = Ax, the 
computer model constructs the spatial frequency meshing. Aft, such that 



ftmax - 2^,, 


(A.2) 


L_sea = As N _sea 


(A.3) 


Ab = 

L sea 


(A.4) 



By invoking the dispersion relation, we convert from the ft domain to the uj domain, using 



UJ = \/2tx gb 



(A.5) 



This process shows how the spatial meshing relates to the sea energy frequency domain 
meshing. As one can see, the fact that N _sea is constrained to be a power of 2 introduces 
some limitations on the mesh size As, since it is required that 



element _length = M As 



(A.6) 



Tow Speed ( [/) 

The tow speed of the antenna is assumed constant with respect to time. Tow speed must be 
greater than zero (in order to generate tension). When 0, the angle between an individual 
propagation direction and the antenna tow angle, is less than 90 degrees (i.e. following seas), 
sea frequency components above cmiega_crit^ 

omega crit = ^ ^ (A. 7) 

are zeroed out to prevent antenna folding. If omega_crit is less than 1.5 rad/sec, then a 
warning is issued that excessive frequency components are being zeroed out, and program 
execution stops. 
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Simulation Time (duration) 

The total length of time simulated by the program. 

Antenna Time Step Increment ( twiner) 

The time step increment used in solving tlie antenna’s governing equation. A typical value 
for an antenna with no stiffness is 0.001 seconds. When stiffness is added, the time increment 
should be reduced to 0.0001 seconds. In all cases, if the program stops executing due to a 
run-time math error, the time increment should be made smaller to allow proper execution of 
the program. 

Snapshot Output Interval (pic^ time) 

This variable allows the user to control the time interval between snapshots, which are recorded 
to file nin_id / /^ _m<xt\ These snapshots record both the sea and antenna elevation data along 
the entire length of the antenna, at a resolution of As. 

Time of First Snapshot (t_ first __ pic) 

The time to begin recording snapshots is user defined. This is useful to ignore transient behavior 
at the start of the simulation, or to save only a representative sample at the end of a lengthy 
simulation. 

Statistics Output Interval (stat^time) 

This variable allows the user to control the time interval between exposure data matrix records, 
which are recorded to file nin_ id/ _edm\ Selection of an excessively small stat_ time results 
in subsequent records containing essentially the same information, and contributes little to the 
statistics of the sample. 

Snapshot Output Option(matlab^output) 

The snapshot output file (nm_ id//’_mat’) is optional, since it does not contain any unique 
statistical information, and requires a relatively large amount of disk space. The primary 
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benefit of examining the snapshot output file is to learn the physical behavior of the antenna in 
a seaway. In addition, the output is useful to detect when the time step size is too small, which 
is characterized by rapid uiiphysical antenna oscillations. rnatlab_output=, TRUE, means that 
snapshot output is enabled, .FALSE, means snapshot output is disabled. 

Simulation Parameters and Results Summary Option {run_info) 

The simulation parameters and results summary file (nin_ z(i//’_info.txt’) is optional, but is 
highly recommended as it contains all of the pertinent inputs and results of the simulation. 
rarz_ m/o=.TRUE. means that summary output is enabled, .FALSE, means summary output 
is disabled. 

Exposure Data Matrix Option {exp_data) 

The exposure data matrix file (ra7i_i(i//’_edm’) is optional; since, it is needed only when 
examining the exposure of individual elements rather than being concerned with the overall 
exposure of the antenna. exp_data=.TKVE. means that the exposure data matrix output is 
enabled, .FALSE, means that the exposure data matrix output is disabled. 

Temporal Wave Data Option {wave_temporal) 

The temporal wave data file (ran_id//’_mat2’) is optional, since it is needed only when an- 
alyzing the power spectrum of the modeled seas in time. It contains the wave elevation at a 
single random point on the antenna in time. wave__ temporal=. TRUE, means that the temporal 
wave data output is enabled, .FALSE, means that the temporal wave data output is disabled. 

Spatial Wave Data Option {wave^ spatial) 

The spatial wave data file (rarz_zc?//’_mat3’) is optional, since it is needed only when analyzing 
the power spectrum of modeled seas along the antenna. It contains the wave elevation at every 
point on the antenna for the time = 0.0 seconds. wave_spatial=,TR\]E, means that the spatial 
wave data output is enabled, .FALSE, means that the spatial wave data output is disabled. 
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Significant Wave Height {sig^wave^ht) 

This variable sets the sea’s significant wave height, a parameter used to describe the sea severity 
in both the Bretschneider and Ochi spectral models. Appendix B correlates sea state to 
significant wave height for a fully developed sea. 

Development {development) 

This variable allows the user to shift the sea spectral energy modal frequency up, for developing 
seas, or down, for decaying seas. This parameter applies only to the Bretschneider spectral 
model. 

Spectral Model {spectral _model) 

Four spectral model subroutines are available: the two parameter Bretschneider spectrum, 
the Ochi “most probable” six parameter spectrum, a square wave spectrum, and a single sine 
wave spectrum. Option’ 1’ means Bretschneider, ’2’ means Ochi, ’3’ means square, and ’4’ means 
sinusoidal. The square wave spectrum and the sine wave spectrum were included for debugging 
and validation of the results. Parameters of the square wave spectrum cannot be changed with 
the input file, but the sine wave amplitude and frequency are computed from the sig_xuave_ht 
and development variables. The user is directed to review the source code for the sine wave 
subroutine to determine how these two input variables control the sine wave parameters. 

Selection of the spectral model has a measurable impact upon the exposure statistics com- 
puted by the code. There is no litmus test that the user can apply in selecting the correct 
spectral model, but the different models are better suited for some specific conditions. The 
assumptions used in deriving the Bretschneider spectrum from North Atlantic sea data were 
deep water, unlimited fetch, uni-directional seas, and no swell. The Ochi spectrum was derived 
using the same assumptions except that it allows for swell from a distant storm. Based on 
those assumptions, one can see that the Bretschneider spectrum is more suited for open ocean 
wind driven seas, and the Ochi spectrum is more suited for open ocean wind and swell driven 
seas. 
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Directional Spreading Function {spreading _^func) 

This allows the user to select the directional spreading function to apply to the power spectra 
to obtain the desired wave field. Two directional spreading functions are available: cos^O and 
cos'^O. Option ’r means and option ’2’ means cos'^O. 

Number of Propagation Directions (NN) 

This allows the user to select the number of discrete unidirectional waves which will combine 
to form the multidirectional sea surface. A larger number of directions will result in more 
realistic sea spectra but at the expense of computation time. It is recommended that the user 
experiment with the number of propagation directions until the simulation results converge. 

Trapezoidal Integration Convergence Factor (epsilon) 

This refers to the convergence factor of the trapezoidal integration that is needed to calculate 
the integral of the directional spreading function. A typical value is le“^. 

Antenna Tow Angle (tow_angle) 

Antenna tow angle with respect to the primary axis of propagation. A tow angle of 0 degrees 
represents following seas, 90 degrees represents beam seas, and 180 degrees represents head 
seas. 

Antenna Diameter (diam_^ antenna) 

This is the outer diameter of the hose-like structure encasing the antenna. Element performance 
(gain) is controlled by whether this outer diameter is submerged; the diameter of the inner 
conducting slot antenna affects the rotational invariance of the element, and the maximum 
theoretical gain of a fully exposed element. 

Water Density (rho_water) 

Water density in which the antenna operates. 
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Drogue Size (drog^area) 

Tlie drogue frontal area (in in^) determines the tension at the trailing end of the antenna, 
according to 



T = ^Pwater^'^ (irog _areaC D _n 



(A.8) 



A typical value is 0.01 m^. 

Normal Drag Coefficient (CD^n) 

This represents tlie drag coefficient for a cylinder in cross flow. A constant value of 1.0 hcis 
been assumed. 

Tangential Drag Coefficient (CD^t) 

This is the skin friction associated with flow along the long a:xis of the antenna. A relatively 
smooth outer surface has been assumed in choosing this value, which has been assumed constant. 
A typical value is 0.0035. 

Antenna Density {rho_antenna) 

The average antenna density, which is assumed to be constant over the entire array. This is 
an approximation in that the antenna has slightly higher density at the actual element, and 
slightly lower density between elements. 

Antenna Lateral Stiffness {flex ^ stiff ^coeff) 

This variable represents the resistance of the antenna to bending, and is given by 



flex_stiff_coeff = El 



(A.9) 



where 



E is the modulus of elasticity 
I is the area moment of inertia. 
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.04 




Figure A-1: Sample Threshold Antenna Exposure of 0.04 
Added Mass Coefficient (C_m) 

The added mass coefficient for a cylinder in cross flow has been set to 1.0. 

Maximum Number of Samples {max^num^ of _ samples) 

(r) 

This refers to the maximum number of ‘‘snapshots” taken for the Matlab^ viewer, or the 
exposure data matrix. max_nu7n_ of _ samples is used to dynamically allocate array storage for 
variables which are sampled periodically in time. 

Threshold Antenna Exposure {reqd^ ant ^exposure) 

As mentioned earlier, the antenna gain remains steady until the element just becomes sub- 
merged. In order to capture this performance, and introduce margin for capillary waves not 
modeled, the parameter reqd_ ant _ exposure has been defined. The reqd_ant_ exposure is the 
height above the water which must be maintained by the antenna to be considered exposed. 
The reqd_ ant _ expo sure is given as a fraction of the antenna diameter. See Figure A-1 for a 
diagram representing an antenna diameter of 1.0 and a reqd_ ant _ exposure of 0.04. 
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Threshold Number of Elements {num_ of __elem__ thresh) 

This variable allows the user to specify the lowest number of exposed antenna elements necessary 
to achieve the desired design antenna array gain. At each stat_ time, the number of exposed 
elements is compared to num _ of _elem_ thresh', the result (a “1” when the threshold is met, 
“0” otherwise) is found in the summary output file, nm_ id//’_info.txt’, along with any other 
variables tracked over the course of the simulation. Currently, the variables tracked over time 
are 



• tow point tension {tow_tension{ )) 

• number of elements exposed {is_exposed{ , )) 

• threshold met {thresh_met: “0” or “1”, as above) 

• counterdetection length {length__ totally _expose(i{ )). 
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A. 2. 2 Sample Input File 



Input File for proteus Simulation Program 

SI units throughout; useful conversions follow 
1 inch = 0.0254 in 
1 foot - 0.3048 m 
1 knot = 0.5144 m/s 



1 5M ay _ 02 


run_ id 


13 


nnf 


60 


num _ of ele 


0.5 


element length 


8 


M 


3.0 


U 


30.0 


duration 


.001 


t _ incr 


0.2 


pic_time 


0.0 


t_ first_ pic 


0.1 


stat_tirae 


FALSE. 


matlab_output 


TRUE. 


run_ info 


FALSE. 


exp_ data 


FALSE. 


wave_tem poral 


FALSE. 


wave_spatial 


3.0 


sig_ wave_ ht 


1.0 


development 


1 


spectral_ model 


1 


spreading^ func 


9 


NN 


le-3 


epsilon 


180.0 


tow _ angle 


0.0762 


diam _ antenna 


1025.9 


rho_ water 


0.01 


drog_area 


1.0 


CD n 


0.0035 


CD^t 


408.0 


rho antenna 


100.0 


flex stiff coeff 


1.0 


C m 


1000 


max _num _of_! 


0.04 


reqd_ ant_ expo 


30 


num _ of_ elem _ 



8 characters long 

# of com ponents. N_sea = 2**nnf 

# of actual slotted antenna elements 
physical slotted antenna length (m) 

# of mesh points/element _ length 
tow speed (use POSITIVE #) (ra/s) 
total simulated run time (sec) 

time step increment used for antenna calculations 
snapshot interval for Matlab viewer (sec) 
time of first snapshot (sec) 
statistics interval (sec) 

TRUE: A file is created that contains wave and antenna data 

TRUE: A file is created that contains a summary of run parameters and results 

TRUE: A file is created that contains the exposure data matrix 

TRUE: A file is created that contains temporal wave elevation at a single point 

TRUE: A file is created that contains spatial wave elevation at time=0 

significant wave height (m) 

modal frequency shift factor (used only for spectral model 1) 

1 = Bretschneider; 2=Ochi; 3 = square; 4 = sine 
directional spreading function: l = cos'2; 2 = cos'‘4 
Number of propagation directions (must be odd number) 

Convergence factor for trapezoidal integration 

Tow angle relative to primary axis of propagation (0-180 degrees)(0 following, 180 head) 

Antenna outer diameter (meters) 

density of seawater (kg/ra''3) 

drogue area normal to fluid flow (m''2) 

normal flow drag coefficient for a cylinder 

antenna/fluid tangential drag coefficient 

antenna density (kg/m "3) 

lateral bending stiffness of antenna (Pa m''4) El 

added mass coefficient for a cylinder 

upper limit on # of samples in time domain 

threshold amount of vertical exposure necess for reception {% of antenna diameter) 
threshold # of elements required to achieve necessary gain 



A. 3 Output Files 

A. 3.1 Snapshot Output File (rtm_ id//’_mat’) 
File Description 
This file provides: 
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Run ID: 24May-01 



Sea Data For Time 0.00 sec 




Figure A-2: Sample Output of Snapshot Output File 

• a time history of sea elevation along the antenna 

• a time history of vertical, longitudinal, and transverse sea velocities along the antenna 

• a time history of vertical antenna position 

(r) 

• antenna parameters necessary for the Matlab^ function file. 

This file is intended to be viewed with the Matlab® function file, proteus_mat. 

Sample Output 

A sample output of the Matlab® function file proteus_mat is shown in Figure A-2 

A.3.2 Simulation Parameters and Results File (rTxn_ 2 d//’_info.txt’) 

File Description 
This file provides: 

• a complete listing of antenna parameters 
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• a complete listing sea spectra parameters 



• a summary of exposure and antenna performance statistics. 



This file is intended to be viewed with any standard text editor. 



Sample Output 



Summary of Parameters and Results 



General Information 
Run ID = 24 May 01 

Parameter and Result Summary File = 24M ay_01 info.txt 
Wave and Antenna Data File = 24May 01_mat 
Start Date/Time = 05/12/2000 07:22:31 
End Date/Time = 05/12/2000 07:41:07 
Elapsed Computational Time = 1115.56 seconds 
Simulated Analysis Time = 30.00 seconds 

Antenna Parameters 



Number of Antenna Elements = 60 

Length of Antenna Element = 0.50 meters 

Total Antenna Length = 30.00 meters 

Antenna Diameter = 0.0762 meters 

Lateral Bending Stiffness = 50.00 N-m"2 

Antenna Density == 408.00 kg/m'3 

Drogue Area Normal to Fluid Flow = 0.0100 m'2 

Normal Flow Drag Coefficient for a Cylinder = 1.000000 

Tangential Drag Coefficient = 0.003500 

Coefficient of Added Mass for a Cylinder = 1.000000 
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Number of Mesh Points Per Element = 8 



Tow Speed = 2.50 m/sec 
Tow Angle = 180.0 degrees 

Time Increment for Antenna Related Calculations - 0.00010 seconds 

Sea Spectra Parameters 

Spectral Model = B retschneider 

Directional Spreading Function = cos' 2 

Development Factor = 1.00 

Significant Wave Height = 3.00 meters 

Minimum Cutoff Frequency = 0.35 radians/sec 

Maximum Cutoff Frequency = 22.21 radians/sec 

Length of Seas Modeled = 512.0 meters 

Number of Spatial Sea Domain Mesh Points = 2" 13 

Number of Propagation Directions = 3 

Seawater Density 1025.90 kg/m'3 

Time Increment for Sea Spectra Related Calculations = 0.02125 seconds 
Number of Antenna Time Increments Between Sea Spectra Calculations 212 

Simulation Results 

Maximum Tow Tension - 116.4 N 
Minimum Tow Tension = 24.3 N 
Stretched Antenna Length = 30.2 meters 
Antenna Strain = .005 

Statistics Sampling Interval = 0.100 seconds 
Number of Samples Observed for Statistics = 301 
Average Fraction of Antenna Exposed = 0.895 

Vertical Antenna Exposure Necessary to Receive 4.0 percent of diameter 
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Number of Antenna Elements Required to be Exposed -30/60 



Time That Number of Required Elements are Exposed 100.0 percent 



Average Detection Length = 1.95 meters 



A. 3. 3 Exposure Data Matrix File {run^idf edm’) 

File Description 
This file provides: 

• the raw exposure data information for each antenna element at each time sampled 

• the spatial position of all elements. 

This allows one to analyze statistics across time or space. This file can be viewed with a 
user supplied Matlab® file or with a spreadsheet application. 

Sample Output 



time 



spatial position / exposure status 



0.00000 30.00000 29.50000 29.00000 28.50000 28.00000 



0.09900 



0 



0.19800 



0 



0.29700 



0 



0.39600 



0 



A. 3. 4 Temporal Wave Data File {run idf mat2’) 



File Description 



This file provides: 
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16 



Power Spectral Density Comparison for Data Run <17May-02> 



14 



12 



10 



CNJ^ 8 



Scaled Output Spectrum 
Input Spectrum 



Area of Input Spectrum = 5.000000 

Area of Unsealed Output Spectrum (Not Shown) = 0.056341 
L/2k Scaling Factor = 86.84 
Area of Scaled Output Spectrum = 4.892339 
Number of data points per segment = 4096 
Number of non-o\«rlapping data segments = 2, 

Total Number of data points = 819^ 




0.6 0.8 1 
frequency (Hz) 



Figure A-3: Sample Output of Temporal Wave Data Analysis 

• a time history of the wave elevation at a single random point on the antenna 

• antenna parameters necessary for the Matlab^ function file. 

This allows one to analyze the power spectrum of the modeled seas in time. It is intended 
to be processed with the FORTRAN executable file, proteus_mat2.exe^ then viewed with the 
Matlab^ function file, proteus_mat2. 

Sample Output 

A sample output of the Matlab® function file proteus_mat2 is shown in Figure A-3. 

A. 3. 5 Spatial Wave Data File (mn_id//’_mat3’) 

File Description 
This file provides: 

• the wave elevation at every point on the antenna at time = 0.0 seconds 

• antenna parameters necessary for the Matlab® function file. 
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16 



Power Spectral Density Comparison for Data Run <l7May-01> 



14 



12 - 



10 - 



Scaled Output Spectrum 
Input Spectrum 



Area of Input Spectmm = 5.000000 

Area of Unsealed Output Spectmm (Not Shown) = 

L/2k Scaling Factor = 64.00 

Area of Scaled Output Spectrum = 4.844453 

Number of data points per segment = 1024 

Number of non-o\er1apping data segments = 8 

Total Number of data points = 8192 

NOTE: Output spectrum data is at time t = 0 



0.075695 



4 

2 

0 

- 0.2 




0.2 0.4 0.6 0.8 

spatial frequency (m*^) 



Figure A-4: Sample Output of Spatial Wave Data Analysis 



This allows one to analyze the power spectrum of the modeled seas in space. It is intended 
to be processed with the FORTRAN executable file, proteus_mat3.exe^ then viewed with the 
Matlab® function file, proteus_ matS. 



Sample Output 






A sample output of the Matlab function file proteus_mat3 is shown in Figure A-4. 
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Appendix B 



Sea State Correlation Data 
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Significant 
Wave Height 
(m) 


Modal 

Frequency 

(rad/sec) 


Wind 

Speed 

(knots) 


Sea 

State 


0.1 


3.96 


4.2 




0.2 


2.80 


5.9 




0.3 


2.29 


7.3 


1 


0.4 


1.98 


8.4 




0.5 


1.77 


9.4 




0.6 


1.62 


10.3 




0.7 


1.50 


11.1 




0.8 


1.40 


11.9 


2 


0.9 


1.32 


12.6 




1.0 


1.25 


13.3 




1.1 


1.20 


13.9 




1.2 


1.14 


14.6 




1.3 


1.10 


15.2 


3 


1.4 


1.06 


15.7 




1.5 


1.02 


16.3 




1.6 


0.99 


16.8 




1.7 


0.96 


17.3 




1.8 


0.93 


17.8 




1.9 


0.91 


18.3 




2.0 


0.89 


18.8 


4 


2.2 


0.85 


19.7 




2.4 


0.81 


20.6 




2.6 


0.78 


21.4 




2.8 


0.75 


22.2 




3.0 


0.72 


23.0 


5 


3.2 


0.70 


23.8 




3.4 


0.68 


24.5 




3.6 


0.66 


25.2 


6 


3.8 


0.64 


25.9 




4.0 


0.63 


26.6 




4.5 


0.59 


28.2 




5.0 


0.56 


29.7 




5.5 


0.53 


31.2 




6.0 


0.51 


32.6 




6.5 


0.49 


33.9 


7 


7.0 


0.47 


35.2 




7.5 


0.46 


36.4 




8.0 


0.44 


37.6 




8.5 


0.43 


38.8 




9.0 


0.42 


39.9 





Table B.l; Sea State Correlation Data 
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